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Abstract 

We analyze the average of weak values over statistical ensembles of pre and post-selected states. The protocol of weak values, 
proposed by Aharonov et al. [3], is the result of a weak measurement conditional on the outcome of a subsequent strong 
(projective) measurement. Weak values can be beyond the range of eigenvalues of the measured observable and, in general, 
can be complex numbers. We show that averaging over ensembles of pre- and post-selected states reduces the weak value 
within the range of eigenvalues of A. We further show that the averaged result expressed in terms of pre- and post-selected 
density matrices, allows us to include the effect of decoherence. 
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1. Introduction 

Non-invasive quantum measurement has been a 
long-standing challenge within the framework of quan- 
tum mechanics, since such a protocol seems to be 
incompatible with the projection postulate [1] in the 
standard formulation of the quantum measurement. 
According to the projection postulate, the measure- 
ment in quantum mechanics is a probabilistic process 
whose outcome is, with a certain probability, one of the 
eigenvalue of the measured observable. Immediately 
following the measurement, the state of the quantum 
system collapses onto an eigenstate of the measured 
operator corresponding to the measured eigenvalue. As 
opposed to the standard picture, non-invasive (weak) 
measurements, weakly disturbs the system, while pro- 
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viding only partial information about the state of the 
latter. A major step in formulating alternatives to the 
standard measurement protocol comes from the two 
state formulation of quantum mechanics, pioneered by 
Aharonov, Bergamnn and Lebowitz [2]. This formal- 
ism was originally introduced to discuss the symmetry 
between past and future during a quantum measure- 
ment, which is explicitly broken by the projection 
postulate. In Ref. [2] it was shown that the result of 
a quantum measurement depends symmetrically on 
both a past state at which the system is prepared 
(pre-selection) and a future state in which the system 
is selected following the measurement (post-selection) . 

Arguably, the most interesting phenomenon in the 
context of two state formalism is the emergence of weak 
values (WV) in measurements between pre- and post- 
selected ensembles [3,4]. The protocol for the weak 
value involves (i) preselection of the system in the 
state Ixo); (ii) meafc measurement of A; (iii) projective 
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(strong) measurement of a second observable B: the 
result of the weak measurement is kept provided the 
measurement of 13 detects the system in a specific, pre- 
selected state \xf)- The weak value resulting from this 
procedure, 

/(^)o = (X/|4|X0>/(X/|X0>, (1) 

can be orders of magnitude larger than the standard 
values [3] , negative [5] (where conventional strong val- 
ues would be positively definite), or even complex. 
Weak values may allow us to explore some fundamen- 
tal aspects of quantum measurement, including ac- 
cess to simultaneous measurement of non-commuting 
variables [6,7], dephasing and phase recovery [8], cor- 
relations between different measurements [9,10], and 
even new horizons in metrology [3] . The access to non- 
commuting variables allows a direct study of entangled 
particles, possibly a direct access to quantum statistics 
of identical particles. 

The realization of a weak value protocol requires a 
high level of control of the dynamics of the quantum 
system at hand and of the process of measurement. 
While some aspects of weak values have been explored 
in optics based experiments [11], the desired high de- 
gree of control is now at hand in quantum solid state 
devices. Indeed, the arena of mesoscopic solid state of- 
fers very rich physics to be studied through weak val- 
ues, as well as the possibility of fine tuning and con- 
trolling the system's parameters through electrostatic 
gates and applied magnetic fields. In fact, the study 
of weak values and their implementation in the con- 
text of solid state system has been initiated very re- 
cently [12,13] and has led to proposal to explore previ- 
ously unexamined aspects of WVs, like their complete 
tomography [14] 

Given that weak values are average values, i.e. are 
obtained by averaging over many repetitions of an ex- 
periment, they will display some noise. Such a noise 
will be two-fold: (i) noise of the weak value keeping 
fixed pre- and post-selected states at a given replica 
of the experiment, of purely quantum mechanical ori- 
gin [15,16]; or (ii) noise due to the effect of replica-to- 
replica fluctuations of the pre- and post-selected states. 
The latter is relevant in the analysis of weak values 
in realistic systems where the pre- and post-selected 
states are obtained by tuning certain system's parame- 
ters; in particular, this is expected to be the case in all 
the reported studies of weak values in optics and solid 



state based systems. Moreover a WV protocol essen- 
tially relies on the quantum coherence of the system at 
hand. Fluctuations of the system's parameters neces- 
sarily affect the measured WVs. 

The present work is devoted to the study of weak 
values averaged over fluctuations in the pre- and post- 
selected ensembles. As we further generalize weak val- 
ues to include pre and post-selected density matrices, 
this sets the ground for studying decoherence effects 
within the WV framework. In the following we present 
the derivation of WV in a simple model of system- 
detector interaction (section 2), and discuss a protocol 
to observe the electron's spin WV in a double quantum 
dot (section 3) as a specific example. We then present 
our results concerning the average of weak values and 
analyze them in the specific example of a double quan- 
tum dot (section 4). We further clarify the relation be- 
tween averages of weak values on one hand, and the 
evaluation of WVs for non-pure states on the other 
hand, and employ it to discuss the manifestation of de- 
coherence effects in weak values. 



2. Weak values 

In an ideal von Neuman measurement, the coupling 
of a system to a detector is described by the Hamilto- 
nian 

H = Hs + HD + Hint, Hint = \g{t)pA , (2) 

where i/s(D) is the Hamiltonian of the system (detec- 
tor), and Hint is the interaction Hamiltonian. Here p 
is the momentum canonically conjugate to the posi- 
tion of the detector's pointer, g, and \g{t) (A <C 1) is a 
time dependent coupling constant. A — "^^i {o-i \ \<^) 
is the measured observable. We assume for simplicity 
that the free Hamiltonians of the system and the de- 
tector vanish and that g{t) — 5{t — to). Before the 
measurement the system is in the state |xo), and the 
detector is in the state |(^o), the latter assumed to be 
a gaussian wave-packet centered at g = go, |(/)o) — 
(7g-(9-9o)V4A2 ^^g^ ^j^g interaction with the detec- 
tor the entangled state of the two is 

It/-) = exp{-iAM} Ixo) |0o) ■ (3) 

A projection onto the state |x/) leaves the detector in 
the state 
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\^) = |0o>-iA(x/|i|xo)/(x/lxo>P |0o) ^ e-'^f^^^o^ |0o) 

(4) 

that corresponds to a shift in the position of the pointer 
proportional to ^e[-^f {A)^„]. Hence the expectation 
value of the coordinate of the pointer (initially equal 
to qo) is given by 

(g) =go-AKe[^(^)^]. (5) 

Under more general conditions (i.e., p, q are not canoni- 
cally conjugated), the imaginary part of the weak value 
may be meaningful too [17] . 

We note that the approximation in Eq. (4) is valid 
if A <C maxij \ai — aj\. This means that the initial de- 
tector's wave function and the shifted one are strongly 
overlapping. In fact, in an ideal strong measurement, 
there is a one-to-one correspondence between the ob- 
served value of the detector's coordinate, q^, and the 
state of the system, \a). Within a weak measurement 
procedure the ranges of values of q that correspond 
to two distinct states of the system, \a) and |a'), are 
described by two probability distribution functions, 
Pa (q) and Pa> (q) respectively. These distributions 
strongly overlap. Hence the measurement of q provides 
only partial information on the state of S. 

3. Weak values of spin in a double quantum dot 

As a reference model we consider a recent proposal 
for observing weak values in a solid state device [12]. 
The device consists of two electron residing on a double 
quantum dot, a setup that has been recently considered 
as a candidate for a qubit [18] . A weak value protocol 
can be realized in this system by means of the (exper- 
imentally demonstrated) coherent control of electrons 
spin and charge distributions in the dots, and the weak 
coupling to quantum point contact to be employed as 
a detector. 

The system (cf. Fig. 1(a)) consists of a gate con- 
fined semiconducting double quantum dot hosting two 
electrons. The charge configuration in the two dots, 
(n_L, Tiij), is controlled by the gate voltages Vl and Vr. 
In particular, by controlling the dimensionless param- 
eter e (X Vr — Vl, the charge configuration is contin- 
uously tuned between (0, 2) and (1, 1). When the two 
electrons are in the same dot (0, 2), the ground state 
is a spin singlet, [^(O, 2)); the highly energetic excited 
triplet states are decoupled. For (1, 1) the degeneracy 




Fig. 1. Schematics of the system and energy levels, (a) 
Scheme of a double dot with nearby quantum point con- 
tacts (QPCs) as charge sensors, (b) Energy levels of low- 
est singlet (red) and triplet (blue) states vs. the detuning 
parameter e. In the (0, 2) charge configuration the anti- 
symmetric nature of the electrons wave function implies a 
singlet ground state. The states \T±) (blue dashed lines), 
with angular momenta itft in the direction of the applied 
magnetic field, are split by the Zeeman energy. The range 
of e in which the effect of nuclear interaction is relevant is 
highlighted by the shadowed part. 

of the triplet states is removed by a magnetic field, 
B — Bz, applied perpendicularly to the sample's plane 
(cf. Fig. 1). Due to the spin-preserving inter-dot tun- 
neling, As(0)/2 (controlled by the gate voltage Vr), 
the ground state of the system is the charge- hybridized 
singlet, \Sg{e)), while the first exited state is |To(l, 1)). 
Theh energy difference, J(e) = As{0)[e + -f 1/4], 
is vanishingly small at e < (cf. Fig. 1(b)). In this case 
the hyperfine interaction between electrons and the nu- 
clear spin [19,20], facilitates transitions between these 
two states. For our purpose, the effect of the nuclear 
spins on the electrons is described by classical magnetic 
fields, Bjv^, Bjvfl, resulting in the Hamiltonian Hn ~ 
gUBiBNR - Bnl) ■ z |To(l, 1)) {3(1, 1)| + H.c. The 
Hamiltonian for the lowest energy singlet and triplet 
states can be then written as 

i? = J(e)|rol,l)(Tol,l| + //iv. (6) 

The state of the system is then controlled by tuning 
the parameter e between e|s and ea- By a variation 
—> fast on the time scale of the nuclear field 
coupling, (ro(l, l)|i?jv|S'(l, 1)), ("fast adiabatic"), 
\S{1, 1)) is mapped into the ground state, [5(0, 2)) and 
[ro(l, 1)) is unchanged. By a "slow adiabatic" varia- 
tion, the ground state of the Hamiltonian Eq. (6), |ti). 
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will be mapped into [^(O, 2)). The described procedure 
is in fact mapping different spin states at e = into 
different charge states at e = ea, and is referred to as 
spin to charge conversion. Notably, by measuring the 
charge of the final state one is effectively measuring 
the singlet, triplet (for "a fast adiabatic variation") 
or the Iti), lit) (for the "slow adiabatic variation") 
component of the initial state at e = es. 

The detectors (Dl and D2) are two quantum point 
contacts (QPCs) located near the dots. They are charge 
sensors [21] suitable for continuous measurements [22]. 
and can be employed both for weak and strong mea- 
surement. D2 is used to perform an effective strong 
measurement of the spin via a spin-to-charge conver- 
sion followed by a strong measurement of charge. By 
contrast, at e ~ es, the charge difference between 
the two spin states is small and a measuring charge 
with Dl corresponds to a weak spin measurement. For- 
mally, the interaction between the double QD and the 
QPC is modeled as Hint = -P(i,i) +^{o,2) -P(o,2), 

where P(„^ is the operator projecting onto the sub- 
space with charge configuration (nL,n/{). de- 
scribes scattering of the electrons in the QPC with 
transmission (refiection) coefficient to (ro)' any incom- 
ing electron in the QPC, |in), evolves to \(t>) = to\t) + 
ro \r), where \t) and |r) are the reflected and trans- 
mitted states for the electron. When the charge con- 
figuration in the double dot is (0, 2), the QPC is de- 
scribed by the Hamiltonian -?/(o,2), corresponding to 
|in), evolves to \ <l)} — t' \t} + r' \r). For e « es, the inter- 
action Hamiltonian can be written as Hint ~ ^^(i.i) + 
J(e)/A,(e)(//(o,2) - i?(i,i)) ® (I - 5V2), where the 
measured observable, A = I - = \Sg(€)) {Sg{e)\, 
is the singlet component of the spin state. Now, if the 
QD is in (1, 1), the Hamiltonian in the QPC is still 
while. If the system is in the \Sg{e)) state, the 
electron in the QPC evolves according to |in) [0') — 
{to + St{e))\t) + {ro + Sr{e))\r) = \cj>) + \A4,). St, Sr c&n 
be tuned to be arbitrarily small in J(e)/As(e). 

The protocol for a weak value is realized a sequence 
of voltage pulses as described in Fig. 2(a). The evo- 
lution of the system in the absence of the detector 
for this protocol has already been realized in experi- 
ment [18]. preselection: Initially, at e = ea, the sys- 
tem is in the ground state, [^(O, 2)). By a fast adia- 
batic variation (cf. Fig. 1) it is evolved into 15(1, 1)) 
(e = es at time t = — r'). This state evolves un- 
der the infiuence of the nuclear interaction until time 
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Fig. 2. (a) A protocol to measure weak values of two 
electron spin: shown are V, V , (voltage bias across 
Dl and D2 respectively) and e. (b) The weak value, 
j(A)q^' = [1 — ^(S'^)q^^]/2, as a function of the parame- 
ters a and 7 for (3 = n. The dark region defines the range of 
parameters for which a positive post-selection is obtained 
with probability ^"(0,2) < 0.5%. The shadowed region (par- 
allel lines) corresponds to the values of the parameters for 
which J (S^)*^' <0. 



t = G, thus preselecting |xo) = cos a \Sg) — isina [To) 
with a = gfisi^NR — Bjv_l) ■ zr'. The evolution of 
the system during the measurement pulse, U{t, 0) = 
\Sg) (Sgl -Fexp(-i/3) [To) (Tol, with /3 = J(e)r, is reab- 
sorbed to define the effective preselection state at t = 



Xo) 



\S,) 



i/3r 



sin a [To) 



(7) 



weak measurement: The interaction of the system and 
the QPC creates an entangled state at time t = r, \ip) = 
W(r,0) |x(j) |0)+iW(r,O) |xo) [Acp) . post-selection: The 
following evolution during the time interval (r, r + r") 
is governed by the nuclear interaction and that from 
T + t" to f is a fast adiabatic variation. At this point 
the state [^(O, 2)) is post-selected by the detector D2. 
The evolution W(f, r) defines the effective post-selected 
state at time t = r, \x'f ) = W"^(f , r) [^(O, 2)) = 
C0S7 \Sg) -I- isin7 [To), where 7 = g^s(Bjvii - Bnl) ■ 
zr" (cf. Fig. 2(b)). 

Measuring the current in the QPC conditional to 
the positive outcome of the post-selection of [^(O, 2)), 
gives, at lowest order in toSt, the weak value of A, 
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/(/>o ~ /o + {2e^V/h)2Re{f{A)[;^h*o5t} , (8) 



ix'Ax'o) 



{x'fMr, 0)\xo) 



(9) 



Here 7o — 2e'^V\tof /h is the current for the (1,1) 
charge configuration. 

With the specific pre- and post-selected states, the 
weak value of U{t, 0)A is then 

f{A)l^^ = cos 7 cos q/(cos 7 cos a — e"*'' sin 7 sin a) . 

(10) 

By tuning the duration of the pulses, one can obtain 
a real WV (e.g. for /3 = tt), which is arbitrarily large 
(e.g. for 7 — Q ^ ^/2) - cf. Fig. 2(b). Improvements 
of the protocol to overcome decoherence effects due to 
the nuclear spins fluctuations are discussed in Ref. [12]. 



4. Weak values averaged over pre- and 
post-selection 

The model discussed in the previous section clearly 
shows that the weak value depends on external pa- 
rameters, namely the angles a and 7 that control the 
preselected and post-selected state respectively. The 
weak value is obtained as an average over many repli- 
cas of the experiment, a natural question is then how 
the weak value is affected by averaging over replica-to- 
replica fiuctuations of the pre- and post-selected states. 
The electron spin weak value case will serve as a refer- 
ence to specifically discuss our more general results. 

Given that weak value with specific pre- and post- 
selected states is already obtained as an average 
over many repetitions of the measurement, averaging 
over pre and post-selected states with distributions 
Qo(lx|o)), Qfilxf)) respectively is obtained through 
either : (i) Repeating numerous times the measurement 
with given \xo) and keeping only the properly 

post-selected states, and determining the weak value; 
then iterating the same procedure with difi'erent pre 
and post-selected states chosen according their distri- 
butions, and determining the averaged weak value; or 
(ii) Choosing \xo) and \xf) according to the distribu- 
tion functions Qo, Qf and performing the weak value 
protocol once; then repeating the experiment many 
times with difi'erent states \xo) and \xf) weighted 
according to the distribution functions Qo, Qf, then 
iterating the procedure and averaging the results. 



Having in mind the protocol discussed in section 3, 
the former case would correspond to averaging Eq. (10) 
over fluctuations of a and 7; admittedly it is experi- 
mentally quite artificial. To be specific, consider a case 
where the post-selected state is fixed and only the pre- 
selected state can fiuctuate, i.e. a fiuctuates. The av- 
erage in (i) corresponds to determining the average 
(weighted with Qo) color of Fig. 2(b) along the line 
7 = cost.. On the other hand, the definition of aver- 
aging according to (ii) is the natural result of an ex- 
periment detecting weak values in the presence of ran- 
dom pre- and post-selection. We discuss the two cases 
separately. Below we employ the following notation for 
probability distributions: P{a,b) is the probability for 
both the events a and b to occur, P{a : b) is the condi- 
tional probability for a to take place given that b did 
take place, P{a, 6 : c) is therefore the probability for 
the occurrence of both a and 6 conditional on the oc- 
currence of c. The following formal relations hold: 

P{a, b) = P(a : b) P{b) = P{b : a) P{a) (11) 
^P(a,6) = P(b) (12) 

b 

P{a : b, c) = P(a, b, c) /P{b, c) = P(a, b : c) /P{b : (3)3) 

The formal expression for the average defined in (i) 
is simply: 

Q/(^)Qo=EE Qo(IXo»Q/(lx/))/(^)o. (14) 

|X0> 

This averaged weak value can lie outside the spectrum 
of the eigenvalues of A. 

In order to properly determine the averaged WV 
defined in (ii), we start by noting that, quite generally, 
the weak value can be written in terms of conditional 
probabilities as: 



Ei aiP{ai, \xf)) _ E, a^P{ai, \Xf) I Ixo)) 



E,-P(lx/),«0 E.n«Mlx/)llxo)) ■ 

(15) 

This equation is easily generalized if the preselected 
state is chosen from a distribution Qo(lxo)), 



x/(^>Qo = 



E,E|xo)a.Qo(|xo))J^(a..lx/)llxo)) 
E.E|,„)Oo(|xo))P(a.,|x/)IIXo» 



(16) 

If the measurement of A is weak, the initial state is not 
disturbed by the measurement and we can write 
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E^("-IX/)IIXo»^Kx/ixo)r (17) 

i 

Y,a,P{ar, \xs) \ |xo» ^ x,(i)xol(X/lxo>|' (18) 

i 

where Eq. (18) is derived by plugging Eq. (17) into 
Eq. (15). Substituting Eq. (17), Eq. (18) into Eq. (16) 
we find 

{Xf\Apo\Xf) 



probability 



(19) 



(X/|Po|x/) 

where po = |xo) Q(lxo)) (xo| is the density matrix for 
the initial state given by the distribution Qo(|xo)). The 
result in Eq. (19) is then the generalization of weak 
values to non-pure states. 

If the final state is chosen according to a distribu- 
tion function Q ( | X/ ) ) i then the average over such a dis- 
tribution function has to be taken separately for the 
numerator and denominator. The previous equation is 
then modified as 





(c) 

















E,E|x,)-P(lx/>,«.)Q/(lx/» 



_ E,E|;,,) E|^o) a.Qo(|xo))P(a»,|x/) : lxo))Q/(lx/)) 

E,E|,,)E|,«>Oo(lxo»P(lx/>,««: lxo»Q/(lx/)9°^ 

In the last inequality we include a distribution function, 
Q(lxo)) for the initial pre-selected state. 

This expression is quite general, it is based on 
probability considerations and can therefore describe 
both strong (projective) and weak measurements. The 
difference between the two is in the specific form of 
the probability distribution functions. In the case of 
a projective measurement of A, P{ai, \xf) : |xo)) = 
Kx/l'^i)Pl('^i|xo)P- In this case, Eq. (20) reduces to 

J2- ai Tr[UaiPa] TrpaiPo] 



Q/(^)Qo =Pf {^)po 



Y:,Tr[Ua,po]Tr[na,po] 



(21) 

where po = |xo> Q(lxo)) {xo| and p/ = |x/) Q(IX/>) (x/l 
define the initial and final density matrix respec- 
tively, and llaj = \ai} {ai\. It appears that f{A}o is 
within the range of eigenvalues of A, and it is in- 
variant under the exchange of |xo) and |x/)- In the 
limiting case of uniform distributions, the final results 

are: Qf^conat{A)Q^ = Tr{Apo}, Q f (A) Qg^const = 
Tr{Apf}, Qj:=const(A)Q(,=const = Tr{A}. 

In case of a weak measurement, the initial 
state is unchanged (to lowest order in the system- 
detector coupling), and we can then assume that 



Fig. 3. Averaged weak values of the electrons spin, (a) 
Sketch of the distribution functions for the parameters con- 
trolling the pre- and post-selected states, a and 7 respec- 
tively. We assume the probabiUty distribution for a to be 
Po(q) = l/(2A)6>(a -I- A)6»(A - a) and the probability dis- 
tribution for 7 to be P/(7) = l/(2r)e(7 -|- r)6»(r - 7). (b) 
Weak value for a preselected ensemble. The plot is the weak 
value as a function of for A = (full line), A = 0.25 
(dotted line), A = 1 (dashed line), A = tt (dash-dotted 
line). The post-selection is on a pure state, i.e. F = and 
70 = 7r/3. (c) Weak values for a post-selected ensemble. 
The plot is the weak value as a function of 70 for F = 
(full line), F = 0.25 (dotted line), F = 1 (dashed line), 
F = TT (dash-dotted line). The preselection is on a pure 
state, i.e. A = and ao = vr/lO (d) Weak values for both 
pre and post-slected ensembles. The plot is the weak value 
as a function of ao for A = (full line), A = 0.25 (dotted 
line), A = 1 (dashed line), A = tt (dash-dotted line). The 
post-selection is on an ensemble described by the probabil- 
ity distribution Pf{'y) with 70 = 7r/3 and F = 0.25. 

E,mX/>,«. : lX0» = P{a.;\xf) : lXo» = 
l(x/IXo)|'. Also j:taiP{\xf),a, : |xo» = Ti^iPiai : 
\Xo) ,\Xf))P(\Xf) ■■ lxo» = /(4>o|(x/IXo)|' = 
{XflMxo) Ix/lxo)- Using these expressions in Eq. (20), 
we find 

Qf {A)qo = Tr{pfApo}/ Tr{pfpo} , (22) 

with the same definition of po, p/ as in Eq. (21). 
Notably, the weak value undergoes a complex con- 
jugation operation under the exchange of Ix),) and 
\Xf), i-e. time reversal symmetry. The general behav- 
ior is reported in Fig. 3; for uniform distributions we 

obtain: Q^=const (^>Qo = Tr{Apo}, Q ^ (A) Q„=conat = 

Tr{Apf} 

1 Qj^— const Qq— const — 

Tr{A}. This demon- 
strates that if we have incoherent distributions for the 
initial or final state, weak values reduce to conven- 
tional (strong) values. 
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As we have already mentioned, the resuh in Eq. (19) 
is the weak value for a given initial density matrix, 
obtained through probabilistic arguments. We may re- 
produce the same Eq. (19) with the measurement pro- 
cedure described in section 2, but replacing the initial 
pure state by a density matrix, po- In this case the to- 
tal density matrix, R of the entangled system-detector 
state after the measurement is given by 

R = po®\4>o) (0o|-jAp[4po®P |(/'o) {(po\-poA®\(j>o) (MP] ■ 

(23) 

Subsequently performing the post-selection, and even- 
tually evaluating the (average) shift of the detector 
wave packet (to first order in A) yields the real part of 
the weak value. Within this approach we can discuss 
the role of decoherence. If we have initially a pure state, 
but the system is affected by fluctuations leading to de- 
coherence during the weak measurement, it will evolve 
to a density matrix. If the fluctuations commute with 
A (i.e. fluctuations is described by a term ^{t)B in the 
Hamihonian with [A,B] = 0) Eq. (19) still holds re- 
placing po with p(r), r is the post-selection time. This 
is 

Xf{A)po = {Xf\Mr)\Xf)/{Xf\p{r)\Xf), (24) 

with p(t) = (W(r,0)poWl"(T-,0))stoc, where W(t, 0) is 
the time evolution operator and the average is intended 
over the the fluctuation of the stochastic parameter. 

If the fluctuations do not commute with the operator 
to be measured, a general expression does not exist, 
but perturbation in the coupling to the fluctuating field 
can be carried out. 



5. Conclusions 

The results presented here clarify the meaning of av- 
eraging weak values over fluctuations of pre and post- 
selected states, which is both a conceptual and exper- 
imentally relevant issue. We have demonstrated that 
averaging, which results from an experimental proce- 
dure, is not simply the average of the weak value over 
the fluctuations of the parameters (as in Eq. (14)), but 
rather a properly deflned weak value for pre and post- 
selected density matrices (cf Eq. (22)). We have also 
shown that the derivation of weak values for a density 
matrix sets the basis for incorporating decoherence ef- 
fects into the weak value protocol. 
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